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(54)'ntle: FILM NOISE REDUCTION BY APPLICATION OF BAYES THEOREM TO POSmVE/NEGATIVE FILM 

(57) Abstract 

At least one frame of a 
strip of photographic film Is ex- 
posed to an image, and, prior to 
being developed, six (or so) cali- 
bration strips in a specially re- 
served frame on the same film 
are individually exposed to six 
predetermined calibration light 
sources producing six progres- 
sively increasing exposure levels 
spanning the exposure latitude of 

the rdm. A film layer develops a csa. ^m^^^ooa ,702 eoo ?w 

positive image in one color and a 
negative image in another color 
using positive and negative dyes, 
respectively. The film is scanned 

to determine the positive and ne- V\\|1^I^ir ^sozb 1 "^"V""' e 300b 

gative dye amounts, Dp and Dq 
at each scanner pixel. The likel- 

iest exposure value E as a func- . • 

tion of Dn and Dp is determined probabilisticaUy from six histograms obtained by scanning the six calibration patches on 
the film, all of the I's thus obtained being stored in a look-up table addressed by the corresponding values of Dp and D^. 
The images recorded on the film are then scanned pixel-by-pixel. tiie Dp. and Dn of each pixel addressing the look-up table 
to obtain the UkcUest exposure value E for that pixel, which is tiien furnished to a printer. By thus combining the posiUve 
and negative images on die film, film noise is reduced in the resulting image to obtain the likeliest exposure value E at each 
pixel with die smallest error possible, with no prior knowledge of die film characteristics. 
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FILM NOISE REDUCTION BY APPLICATION OF 
BAYES THEOREM TO POSITIVE/NEGATIVE FILM 

5 TECHNICAL FIELD 

This invention is related to processes for 
improving images obtained from photographic film and 
in particular to reducing film noise in the 
photographic images so obtained. 

10 BACKGROUND OF THE INVENTION 

Main Problem 

Noise attributable to innate 
characteristics of photographic film limits the 
quality of images obtainable from photographic 

15 film. The major source of such noise is random 
variations in the density of silver halide grains 
across the surface of a film emulsion layer. If the 
film is divided into many picture elements 
("pixels") which are sampled individually / each 

20 pixel will have a different number of silver halide 
grains. Such random variations in the number of 
silver halide grains in each "pixel" on the film 
causes different pixels to respond differently to 
the same incident. light . Upon chemical development/ 

25 pixels with more grains develop higher dye 
= densities. These pixel-by-pixel variations in 

density appear as granularity. This can be thought 
of as random pixel-to-pixel variations in 
sensitivity — pixels with more grains are more 

30 sensitive, having a higher response to the same 
incident light. 

However, by independently developing and 
digitally scanning both latent and non-latent silver 
halide grains in each pixel, the total number of 

35 grains in each pixel can be determined. Using this 
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inforrnation, the random variations in the number of 
grains per pixel on film may be compensated by 
simply normalizing or dividing the number of latent 
(or non-latent) grains counted in a given pixel by 
5 the total number of latent and non-latent grains in 
that pixel • 

Typically, using a single dye (positive or 
negative), scanning the film yields either the 
number of non-latent grains or the number of latent 

10 grains in each pixel, but not both numbers. In 
order to obtain both numbers, two different color 
dye are used, one dye responsive only to the 
non-latent grains and the other dye responsive only 
to the latent grains. By developing the same film 

15 with both dyes, and by scanning the film with two 
different color filters tuned individually to each 
of the two dyes, the scanner can separately count 
the number of latent and non-latent crystals in the 
film. Such a process is described in detail by 

20 Bird, "Normal Development, Reversal Development, and 
Composite Processing: A New Method for Gaining a 
Simultaneous Improvement in Latitude and Detective. 
Quantum Efficiency in Silver Halide Films," 
. Photographic Science and Engineering . Vol. 22, No. 

25-6, pages 328-335, November /Dec ember, 1978 

(hereinafter referred to as "Bird"). The Bird 
reference proposes an algorithm for combining the 
densities of D^ and D^ by individually weighting 
them and then computing the exposure from the 

30 combination. The Bird algorithm optimizes the 

detective quantum efficiency (the ratio between the 
input noise and the output noise of the process) 
requiring certain assumptions in film 
characteristics, the main assumption being an 

35 idealized monodisperse emulsion in the film. This 
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assumption means that the silver halide grains are 
of identical size throughout the emulsion so that 
each grain: (1) presents the same photon absorption 
cross-section "area"/ (2) requires the same number 
5 of photons in order to become latent, (3) lies at 
the same height as all other grains in the emulsion 
so as not to obscure one another, and (4) produces 
the sjame amount of dye. 
Secondary Problem 

10 The foregoing assumption clearly cannot be 

relied upon in processing real film. Therefore, the 
process, as described in the Bird reference above, 
introduces errors or noise into the resulting 
reconstructed image • This is because the process 

15 does not take into account secondary sources of 
noise in real film, specifically the 
non-monodisperse nature of the real film emulsion. 
The grains in such a real emulsion do not have the 
same photon absorption cross-section, do not require 

20 the same number of photons for becoming latent, do 
not lie at the same height in the film layer, and 
produce different amounts of dye. Thus, known 
processes, such as that described in the Bird 
reference, for reducing film noise introduce 

25 additional sources of error attributable to 
secondary noise sources in the film. 

PROBLEM TO BE SOLVED 
What is needed is a process for reducing 
film noise attributable to random variations in the 

30 density of active grains across a film emulsion 
without requiring prior assumption of film 
characteristics, and which is therefore less 
affected by any of the sources of film noise 
discussed above- In the present invention, it is 

35 recognized that the prior approach to solving this 
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class of problems was unworkable. This approach was 
based upon trying to understand the sources of noise 
in film, to mathematically model them, and then by 
some mathematical technique invert the process to 
5 remove the noise. The Bird reference is an example 
of this approach in which only the primary source of 
noise was considered while the secondary sources 
were ignored. 

In the present invention, a different 

10 approach is taken, namely to eschew any 

understanding of the sources of noise in film or any 
attempt to mathematically model them, and instead 
process the data from the positive and negative dye 
images in such a manner as to automatically take 

15 into account all sources of noise without having to 
understand them. The advantage of such an approach 
is that it does not depend upon the simplifying 
assumption of the algorithm proposed in the Bird 
reference (namely the assumption of a monodispersive 

20 emulsion), and therefore can be used with a real 
film. 

SUMMARY OF THE INVENTION 
A user exposes the frames on a roll of film 
to a number of scenes, but the last frame on the 

25 V. roll of film is reserved for calibration purposes. 
Before the roll of film is developed, the last frame 
is divided into six patches and the six patches are 
exposed to six different predetermined exposure 
levels comprising a set of six exposure levels 

30 spanning the estimated film latitude. (The number 
of calibration patches is arbitrary.) If the film 
comprises more than one emulsion layer (such as 
color film), then one (or more) emulsion layer is 
developed using positive and negative dyes of 

35 different colors responsive to the non-latent and 
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latent grains in the film, respectively. A scanner, 
such as a CCD imager/ scans each of the six 
calibration patches twice, on the first pass through 
a filter tunefl to the color of the positive dye and 
5 on the second pass through a filter tuned to the 
color of the negative dye. 

During the first pass, the scanner measures 
the dye amount from the non-latent grains in each 
pixel in each calibration patch, and during the 

10 second pass measures the dye amount from the latent 
grains in each pixel in each calibration patch. 
From the dye amount measurements thus produced by 
the scanner, six histograms are constructed, one for 
each of the six calibration patches. Each histogram 

15 represents the distribution of the number of pixels 
as a function of the dye amount per pixel of 
positive and negative dye sites sensed by the 
scanner. 

Each of the histograms is treated as a 
20 graph of the probability that, given a particular 
exposure level E, any pixel will develop a certain 
positive dye amount Dp and a certain negative dye 

amount D . The six histograms, which will now be 

n ' 

referred to as probability distributions, are stored 
.25 memory. 

Each user image on the film is now 
scanned. The most probable exposure value for each 
pixel in a user image is computed from the six 
probability distributions previously obtained. The, 

30 following is one way for performing the computation 
in accordance with the invention: The positive and 
negative dye amounts and of a particular 
user image pixel sensed by the scanner are 
observed. Then, the six probability distributions 

35 are each evaluated at the observed values of 
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and Dp to yield six corresponding probabilities* 
Each of the six probability values thus obtained is 
then multiplied by the exposure value of the 
corresponding one of the six calibration patches to 
5 obtain a weighted exposure value. The six weighted 
exposure values thus obtained are added together and 
the resulting sum is divided by the sum of all six 
probability values to obtain the most probable 
exposure value of the image pixel in question. In 

10 proper statistical language, this is not precisely 
the "most probable" exposure, but rather is the 
"mean" exposure, which often is more desirable. The 
foregoing is repeated for each pixel in a frame on 
the film containing the exposed user image. The 

15 likeliest exposure values thus computed are fed to a 
printer which generates the user image. 

The foregoing process is a simplification 
of the general process of this invention. 

In the preferred embodiment of the 

20 invention, a more sophisticated algorithm is used 
for accurately interpolating a continuum of exposure 
values among the six calibration exposures, which 
takes into account the different probability 
distributions involved. Furthermore, an a priori 

25 probability that a pixel will receive a certain 
exposure E is determined and accounted for in the 
preferred algorithm. 

GAUSSIAN INTERPOLATION 
For computation of exposure values lying 

30 between successive ones of the six calibration 

exposure levels, the preferred embodiment of the 

invention includes a unique interpolation process 

which is as follows: 

The mean values of D and D , their 

n p 

35 variances and their covariances are computed for 
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each of the six probability distributions derived 
from the six calibration patches. If the scanner is 
an eight-bit CCD Scanner (for example), a likely 
choice is to divide the exposure scale defined by 
5 the six calibration exposures into 256 exposures. 
The mean/ the variance and the covariance of 
and of Dp are each plotted against the six 
calibration exposure values and interpolated at 256 
points along the exposure scale. Preferably, a 

10 well-known method of interpolation using cubic 

splines is employed. This produces 256 interpolated 
values of the mean, the variance and the 
covariance. From these, 256 new probability 
distributions are constructed from each of the 256 

15 sets of mean, variance and covariance values by 
simply assuming that the new probability 
distribution at each one of the 256 interpolated 
points is multivariate-Gaussian. Each of the 
Gaussian probability distributions is, in accordance 

20 with well-known principles, uniquely defined by the 
corresponding mean, variance and covariance at the 
corresponding one of the 256 interpolated exposure 
values. 

Next, the image processing algorithm 
25 described in the previous sub-section above is 

performed but, instead of using the six probability 
distributions obtained from the six calibration 
patches, the 256 new probability distributions are 
used. The likeliest exposure value for a given user 
30 pixel is computed as a weighted sum over all 256 
probability distributions, resulting in a 
"likeliest" exposure value. 
Exposure Probability Distribution 

In a further refinement of the preferred 
35 embodiment of the algorithm, each of the probability 
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distributions, corresponding to a given exposure 
level, is first multiplied by the probability of 
receiving that exposure value before it is used in 
the algorithm. The probability of receiving a given 
5 exposure level is computed a priori from a large 
collection of representative photographs or images. 
For example, a large collection of consumer 
photographs may be scanned if the user image is to 
be processed for consumer type images. On the other 

10 hand, in order to process astronomical images, a 

large collection of astronomical photographs may be 
scanned to compute the required exposure 
probabilities. 

SCRATCH DETECTION AND SUPPRESSION 

15 The algorithm may also be used to locate 

defective pixels, such as from scratches or dust. 
Whenever a particular pixel is detected with joint 
positive and negative dye amounts D^ and D^ 
yielding a very low probability when substituted 

20 into each of the probability distributions (such as 
those constructed from the calibration patches), 
that pixel is designated as a defective pixel. 

A scratch suppression algorithm may then be. 
employed. In one such scratch suppression 
. ; 25* algorithm, ; for each pixel recognized as being 

defective, the suppression algorithm defines a small 
neighborhood of pixels around the defective pixel. 
It searches that neighborhood for any non-defective 
pixels, and keeps enlarging the neighborhood until 

30* it finds at least one non-defective pixel in the 
neighborhood. It then computes the median of all 
the non-defective pixels in the neighborhood, and 
replaces the original defective pixel by this median, 
ADVANTAGE OF THE INVENTION OVER THE PRIOR ART 

35 The probability distributions obtained from 
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the calibration patches characterize all noise 
sources present in that film, as well as in the 
scanner. These noise sources include random 
variations in grain density, random variations in 
5 grain size, random variations in photon absorption 
cross-section per grain, random variations in number 
of photons required to transform a grain to 
non-latency from latency, and variations in grain 
height in the emulsion causing some grains to 

10 obscure other grains, and even effects occurring 

only after development such as dye clouds obscuring 
other dye clouds and variations in the dye per 
grain. To the extent that the film response is 
statistically uniform everywhere on the strip of 
.15 film, the probability distributions obtained from 
the calibration patches accurately reflect all of 
the foregoing noise sources. Thus, the image 
processing algorithm of the invention requires 
absolutely no knowledge whatsoever of the film 

20 characteristics prior to execution of the 
algorithm, (The only exception is that the 
progression of exposure values of the plural 
calibration patches should span the exposure 
latitude of the film.) Thus, the algorithm makes no 

25 assumptions about grain size distribution, 

variations in the number of grains per pixel, fog 
levels, stain levels, non-linearity in the 
development process, photon thresholds of the grains 
or many other sources of noise. It does not depend 

30 upon any mathematical model to explain the behavior 
of the film, but instead relies on the calibration 
patches to provide an empirical "understanding- of 
the film. 

It can be shown mathematically (based on 
35 Bayes' Theorem) that the algorithm of the present 
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invention provides the estimate with the minimum 
mean squared error for the exposure level o£ each 
pixel, and any other algorithm will result in a 
larger mean squared error* It is in this sense that 
5 the algorithm o£ the invention is optimal. 

DESCRIPTION OF THE DRAWINGS 
The invention may be understood by 
reference to the accompanying drawings of which: 

Fig. 1 is a block diagram illustrating the 
10 joint positive and negative development process of a 
film bearing a series of calibration patches; 

Fig. 2 is a diagram illustrating the 
generation of an individual probability distribution 
from each one of the calibration patches on the film; 
15 Fig. 3 is a graph illustrating the 

cross-sections of the probability distributions 
obtained from all the calibration patches on the 
, film; 

Fig. 4 is a block diagram illustrating the 
20 image processing algorithm of the invention using 
the probability distributions of Fig. 2; 

Fig. 5 illustrates a refinement of the 
algorithm of Fig. 4 in which the several probability 
distributions of -Fig. 2 are interpolated to generate 
'r 25 many probability distributions; 

Fig. 6 illustrates the type of results 
obtained from the interpolation of Fig. 5; 

Fig. 7 is a block diagram illustrating the 
preferred embodiment of the invention; 
30 Figs. 8a-8d illustrate cross-sections at 

four selected amplitudes of the probability 
distributions obtained from histogramming the 
calibration patches on a real film; 

Figs. 9a-9d illustrate the actual results 
35 obtained by interpolating between the probability 
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distributions of Figs. 8a-8d; and 

Figs. lOa-lOd illustrate the interpolated 
probability distributions of Figs, 9a-9d at the same 
exposures as Figs. 8a-8d to demonstrate that the 
5 interpolated data resembles the original data, 
DETAILED DESCRIPTION OF THE INVENTION 
Film Calibration 

Referring to Fig. 1, a roll of film 100 may- 
have, for example, 36 frames. In the present 

10 invention, the user is allowed to take pictures on 
the first 35 frames, the 36th frame being reserved 
for use by the film development lab and is therefore 
not exposed by the user. (Alternatively, it may be 
exposed during manufacture, or by a special device 

15 within the camera.) When the film 100 is received 
by the film development lab, frame 36 is divided 
into a number of patches (for example, six patches), 
each patch being exposed to a different 
predetermined exposure level from a calibrated light 

20 source 105. One way of doing this is to simply lay 
a six exposure level gray scale mask 110 over frame 
36 and then activate the calibrated light source 
105. 

Joint Positive and Negative Image Development 
25 The film 100 is then developed using a 

special joint positive-negative image development 
process. If the film 100 is, for example, color 
film, it essentially comprises three different color 
layers. In one realization of this invention, only 
30 oiie film layer (usually the top, green layer) is 
subjected to the joint positive-negative image 
development process, the remaining two layers being 
developed using ordinary processing techniques. The 
joint positive-negative image development process of 
35 the invention is illustrated in Fig. 1 and is as 
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follows : 

The top film layer is first developed by a 
partial latent grain developer 120 which develops 
only a fraction of the latent grains in the top 
5 layer of the film 100. Latent grains, as the terra 
is employed in this specification, are those silver 
halide crystals which have absorbed enough photons 
during the exposure of the film to experience a 
chemical change. The partial grain developer 120 

10 creates a dye of one color (for example, cyan), 
which colors the film at each latent grain site 
developed by the developer 120. The result is a 
negative dye density which varies across the 
top layer of the film 100 so as to define a negative 

15 cyan color image. Thereafter, a black and white 
latent grain developer 122, creating no color dyes, 
develops the remaining latent crystals which were 
not developed by the partial latent grain developer 
120. Thus the latent grains have been developed to 

20 yield a negative dye density defining a cyan 
negative image. 

The film is then processed with a 
non-latent grain developer 124 which develops all of 
the non-latent grains in the top layer of the film 

25 100 creating a different colored positive dye (for 
example a magenta dye) so that all of the non-latent 
grain sites in the top layer of the film 100 are 
marked by a magenta dye spot. Non-latent grains, as 
the term is employed in this specification, are 

30 those silver halide crystals in the top layer of the 
film 100 which have not absorbed a sufficient number 
of photons during the exposure of the film to 
experience a chemical change to make them responsive 
to the previous developers 120 and 122. The result 

35 is a positive dye (magenta) density D which 
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varies across the top layer of the film 100 to 

define a magenta colored image. 

Thus, the top layer of the film 100 bears 

two images of different colors, a negative cyan 
5 image and a positive magenta image. The result is 

illustrated in the developed film strip 100' in Fig. 

7. Each of the six calibration patches in frame 36 

also has a cyan negative image and a magenta 

positive image of generally uniform density across 
10 the patch. Generally, as the cyan density increases 

from patch to patch, the magenta density decreases. 
The problem now is how to overcome the 

effect of film granularity which degrades the 

quality of the user images in frames 1-35. 
15 Computing the Conditional Probability 

Functions of Joint Dye Densities for 

Each Exposure Level 

Each of the six calibration patches in 

frame 36 is now used to construct the probability 

20 function of the joint positive and negative dye 

densities D , D for each of the six exposure 
n p 

levels. This is done as follows: An electronic 
image scanner 200 such as a CCD imager scans each of 
the six calibration patches in frame 36, one at a 

25 time. The scanner 200 views each calibration patch 
as simply a plurality of pixels, the number of which 
depends upon the resolution of the scanner 200. The 
scanner 200 first views the calibration patch in 
question through a green filter (complementary to 

30 magenta) so as to detect the positive dye density 

' D at each pixel. Then the scanner 200 views the 
P 

calibration patch through a red filter 
(complementary to cyan) so as to detect the negative 
dye density D^ at each pixel. The scanner 200 
35 observes each of the positive and negative dye 
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densities D„ as the brightness of a 
p n 

particular pixel when viewed through a green or red 

filter respectively. The data thus obtained from 

the scanner for each calibration patch (exposure 

5 level) is combined as a histogram for that exposure 

level, defining the number of pixels sensed by the 

scanner 200 having a particular joint density value 

Dp, D^. The results are plotted in Fig. 2, in 

which the data created by the scanner 200 from each 

10 of the six calibration patches of frame 36 yields a 

three-dimensional "bell" curve, the two horizontal 

axes corresponding to values of D and D 

n p 

respectively, and the vertical axis corresponding to 
the number of pixels at each joint value D^, 
15 In this invention, each of these histograms 

is viewed as the conditional probability 
distribution for finding a joint dye density value 
D , D given a particular exposure value E. . 

It ^ X • 

In this example, there were six calibration patches 

20 so there are six conditional joint probability 
distribution functions, as illustrated. 

^ A cross-section 3-3 of each of the six 

joint probability distribution functions of Fig/ 2 
may be taken at a- chosen amplitude, say 10% of the 

,25r maximum number of pixels counted by the scanner 
200, Each of these cross-sections may then be 
superimposed as illustrated in Fig. 3, in order to 
provide an understanding as to how the joint 
probability distributions differ among the six 

30/ exposure values of the calibration patches of frame 
36 in Fig. 2. 

The scanner 200 may be thought of as a CCD 
imager characterized by a square array of 
photon-collection sites, comprising n rows and n 

35 columns of such sites. Each site defines one pixel. 
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2 

SO that there are n pixels in the imager. In one 

embodiment of the invention, the scanner takes one 

-snapshot" of each calibration patch through the red 

filter and then one shot of each calibration patch 

5 through the green filter to generate the data 

required. If the signal created at each of the 

photon-collection sites in the scanner 200 is 

converted by an analog-to-digital converter to an 

eight-bit byte, then each pixel has a dynamic range 

10 from 0 to 255 units. In this case the units specify 

the minimum and maximum detectable values of the 

negative and positive dye densities D^, D^. 

Thus, each of the six pairs of horizontal axes 

labeled D and D of Fig. 2 extend from 
n p ' 

15 normalized dye density values 0 through 255, The 

vertical axis represents the percentage of the total 

number of pixels. In this example, the total number 
2 

of pixels is n . Typically, m CCD imagers 

2 

suitable for the present application, n is on the 

20 order of half a million pixels. 

The data representing each of the six joint 
probability distributions illustrated in Fig. 2 is 
stored in memory. In this application the data 
characterizing each of the six joint probability 

25 distributions illustrated in Fig. 2 may be thought 

of as an XY addressable memory, in which the X axis 

corresponds to the 256 possible values of D^, the 

y axis corresponds to the 256 possible values of 

D and the data stored in each of the X, Y 
P 

3.0 locations in the memory are the corresponding pixel 
counts from the scanner 200. Such an X, Y array of 
data shall be denoted in this application P<Dp/ 

I E^) , using conventional notation from 
statistical mathematics, P is the conditional 

35 probability value that, given a particular exposure 
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value E^, a pixel selected at random will be found 
to have a particular joint positive and negative dye 
value D^, Each joint probability 
distribution function illustrated in Fig, 2 is 
5 labeled using this same notation. 

Computation of the Probability of Having 

Received a Particular Exposure Value Given a 

Certain Joint Dye Density Value D , D 

n p 

10 The six joint probability distributions 

P(Dp,D^|E^) generated for the six calibration 

exposures may be used in combination to infer the 

related conditional probability P(E. |D , D ). 

1 ' p n 

In accordance with conventional notation, the latter 
15 is the probability that, given an observed joint 
positive and negative dye density value D^, D , 
a particular exposure value E^^ was received by the 
pixel in question. Given that a particular joint 
dye density value D^, D^ is observed at a given 
20 pixel, the question is: what is the probability 

that that particular pixel received a certain 
exposure value E^^? This quantity is 

The answer to this question is obtained 
25; using the well-known Bayes theorem of statistical' 
mathematics. One additional piece of information is 
needed, the a priori probability for a given 
exposure E^^, designated P{E^) . This is the 
probability that the pixel has received exposure 
30 E^, based on no knowledge of D^, D^. This may 
be obtained by histogramming many typical 
photographs. Then according to Bayes theorem, the 
probability that one of the six exposures E^, ... , 
E was received by the pixel is: 
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6 

P(Ej|Dp, Djj) = P(Dp, DniEj) P<Ej)/i: P(Dp, Dn|Ei) (1) 

i«l 

By interpolating 256 exposures among the measured 
six exposures (Fig. 5), exposure may be 
determined with greater resolution. 
5 This result may be used to enhance the user 

images stored in frames 1-35 of the film 100/ as 
follows : 

Computation of the Most Probable 

10 Exposure Value from a Given Joint Dye Density 

Value D„, D^ 
n p 

The most probable (or more accurately the 
"mean") exposure value may be computed for each 
joint dye density value D^, D^ observable by the 
15 scanner 200, and this value shall be labeled herein 

E (Dp/ ' simply E. 

For a particular pixel, E is the normalized 
sura of the six calibration exposure values E^ (for 

20 all i from 1 through 6) individually weighted 
(multiplied) by the- probability amplitudes 
P(E^|D^,Dp) (as defined above in equation 1) 
evaluated at the Observed joint dye density value 
Dp, D^ of the particular pixel. Using the 

25 definition of equation 1, this weighted sum is: 

_ m 
(2) E (Dp, Dn) « E Ei P(Dp, Dn | Ei) P(Ei)/ 

i«l 

m 

Z P(Dp, Dn I Ei) P(Ei). 
30 i = l 

The evaluation of the six probability 

2 

distribution functions at each of the 256 
possible joint density values Dp, D^ (assuming 
the scanner 200 has an eight-bit resolution) 



35 
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requires 256^ different computations using 

equation 2. The result is 256 values of E 

associated with corresponding ones of the 256^ 

possible joint dye density values D , D . As 

p n 

5 illustrated in Fig. 4, these results are 

stored in a look-up table in which a most probable 

(mean) exposure value E is stored in a location 
addressed by the corresponding joint dye density 
value D^. The look-up table 300 is the key 
10 in processing each of the user images in frames 1-35 
of the film 100 through the scanner 200, For better 
precision, the six exposures E^ are replaced by 
the 256 interpolated exposures. 

Joint Positive-Negative Dye Image Enhancement of 

15 User Images 

Now the scanner 200 scans a given one of • 
the user images stored in one of the frames 1-35 of 
the film 100. Each pixel in a given user image is 
determined by the scanner 200 to have a certain 

20 positive dye density D^ and a certain negative dye 
density D^, each of these dye densities lying 
somewhere between a value of 0 and 255, assuming the 
scanner has eight bits of resolution. For each 
pixel, the scanner 200 transmits the observed joint 

25 dye density value D^, D^^ to the address input 

300a of the look-up table 300. The look-up table 
300 responds (in the manner of a programmable read 
only memory) by transmitting the corresponding 

likeliest exposure value E at its data output 300b. 
A succession of likeliest exposure values defining 
an enhanced version of the user image is thus 
transmitted in an order determined by the scanner 
200 to an output device such as a printer. 

REDUCTION OF NOISE 

35 
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It is clear that obtaining independent 
images from the latent and non-latent grains in the 
film provides more information from a given image 
exposure than development of only one or the other 
5 of the positive and negative images. With this 
additional information, it is clear that the noise 
in the resulting image should be reduced. The 
question is how to combine the information. from the 
positive and negative images without introducing new 

10 sources of noise or distortion. It can be proven 
mathematically that the algorithm described above 
computes a probable exposure value for each scanned 
pixel in the film with the minimum mean squared 
error theoretically possible, and that any other 

15 method will result in a larger mean squared error in 
the computed exposure value. Equations 1 and 2 are 
the direct result of Bayes theorem in statistics, 
and may be simplified as follows: 



20 



(3) P (E|Dp,D^) \ P(I>p'I^nl^> 



1 



P' 



(4). . ^ = f P<°p'°nl^i> 

I ^<'^p'°nlEi) 

25 The simplification of equations 3 and 4 from 



equations 1 and 2 is obtained by assuming P(E) 
(which is the a,, pri ori probability that a given 
exposure level will be received) is the same for all 
exposure values — and is therefore conveniently 

30 defined to be unity. This results in only a small 
degradation in the performance of the algorithm. 

In a more complex embodiment of the 
invention, P(E) is not uniform and is inferred from 
a large collection of photographs from consumer 

35 images {if consumer film is to be processed) or 
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astronomy images (i£ astronomical photographs are to 
be processed), etc. 

GAUSSIAN INTERPOLATION OF THE 
JOINT PROBABILITY DISTRIBUTIONS 
5 From equation 2, it is apparent that in 

most cases the likeliest exposure value E will lie 
somewhere between various ones of the six 
calibration exposure values (used to generate the 
calibration patches of frame 36 of the film 100), 

10 and therefore should be properly interpolated. 

Experimental results obtained using this invention 
indicate that the joint probability distribution 
functions such as those illustrated in Fig. 2 
(obtained from the calibration patches) are 

15 approximately Gaussian in nature. Therefore, for a 
more accurate computation of likely exposure values 

E which lie between adjacent ones of the six 
calibration exposure values, a large number of 
"interpolated" probability distribution functions, 

20 lying (in the D^, plane) between the six 

probability distribution functions illustrated in 
Fig. 2 are generated using an interpolation method 
based upon the Gaussian nature of the probability 
/ distribution functions of Fig 2. 

25 The interpolation process is illustrated in 

Fig. 5 while the results are partially illustrated 
in a cross-sectional view in Fig. 6 corresponding to 
that of Fig. 3. In Fig. 6, the six probability 
distributions generated from the six calibration 

30 patches are illustrated in cross-section in solid 
line. An interpolated probability distribution 
lying between each adjacent pair of solid line 
probability distributions is illustrated in dashed 
line in Fig. 6, the interpolated (dashed line) 



35 
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probability distribution having been obtained using 
the interpolation process of Fig. 5, 

In the preferred embodiment of the 
invention, 256 interpolated probability distribution 
5 functions (not shown in Fig. 6) are generated 

corresponding to 256 exposure values lying in the 
range of the minimum and maximum calibration 
exposure values and Eg. Preferably, E^ and 
Eg substantially span the exposure latitude of the 
10 film 100. 

As illustrated in Fig. 5, the following 
process is performed for each of the six probability 
distribution functions illustrated in Fig. 2. 
First, the mean value of the negative dye density, 

'^^ D„(E.) at the i^*^ one of the six calibration 
exposures is computed from the i probability 
density function in accordance with equation (a) 
illustrated in Fig. 5. Six such mean values 

D^(E^) are thus obtained for each of the six 
2q values from of i from 1 through 6. When these six 
values are plotted, the result is the graph 

illustrated in Fig. 5 of vs. E. The same 

n 

process is performed using equation (b) to compute 

25 the six mean values Dp(Ej^)/ which generate six 

data points illustrated in the corresponding graph 

of D vs. E in Fig. 5. The variance S (E.) 
p nn 1 

is computed in accordance with equation (c) in Fig. 

5 for each of the six probability distribution 

30 functions of Fig, 2, resulting in the corresponding 

graph of S^^ vs. E in Fig. 5. In accordance with 
nn 

equation (d) , the variance of the positive dye 
density values ^pp^^i^ similarly is computed, 
which results in the graph of S^^ vs. E in Fig. 

35 
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5. Finally, the covariance S^^ of the positive 

and negative dye density values is computed in 

accordance with equation (e) of Fig. 5 for each of 

the six probability density functions of Fig. 2, 

5 resulting in the corresponding graph in Fig, 5 of 

S vs. E. 

pn 

If, for example, it is desired to compute 

the likeliest exposure value E with eight-bit 
resolution, then each of the graphs of Fig, 5 is 
interpolated so that 256 data points lie in the 
range from iai through i«6. The preferred method 
for interpolation is interpolation using quasi-cubic 
Hermite splines. The method is well-known and 
commercial software is available for performing the 

15 interpolation by cubic splines. The goal in 

interpolating between the six data points in each of 
the graphs above is to achieve the smoothest 
possible curve which most truly represents the six 
data points. The horizontal axes of each of the 

20 five graphs of Fig. 5 may now be thought of as 
having been divided into 255 equal segments 
separated by 256 points, with an index j which runs 
from 1 to 256, 256 points lying between the points 
designated i«l and i=6. For each value of j, the 

25 corresponding values of D (E. ) / D (E. ) / 

off" the interpolated curves and a Gaussian 
probability distribution is then constructed from 
these five particular values in accordance with 
30 equation (f) of Fig. 5. This process is repeated 
for each value of j, so that 256 probability 
distributions are computed in accordance with 
equation (f ) , These are the interpolated 
probability distributions. It is the interpolated 
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probability distributions (from equation f of Fig, 
5) which are now substituted into equation 2 in 
place of the six original probability distributions 
(computed directly from the six calibration patches) 
5 to increase the exposure value resolution of the 
algorithm of Fig. 4. With this change, the 
summations in equation 2 (and in Fig. 4) are over 
256 points (namely from i=l to i=256) instead of 
only six points as indicated in equation 2. 
10 Equation f in Fig. 5 is the well-known 

definition of the multivariate Gaussian probability 

distribution for two variables, in this case D 

-1 " 
and D . The superscript in equation f 

. T 
indicates the inverse of a matrix, while the 

15 superscript indicates the transpose of a vector, in 
accordance with standard notation, 

SCRATCH SUPPRESSION 
Many film defects, such as scratches, 
pinholes and dust, result in unreasonable values of 

20 (Dp, D^) • More precisely, "unreasonable" means 
that the defect produces a joint density (D^, 
D^) for which the probability P(Dp, Dj^lE) is 
very small for all possible values of E. Thus, by 
slightly modifying the algorithm to watch for 

25 unreasonable values of (D^, D^) , the algorithm 
can recognize many film defects. 

A scratch detection algorithm and a scratch 
suppression algorithm have been added (here the word 
"scratch" is meant to be a generic word for the 

30 defects that can be detected) , which recognize 

scratches as described above and suppress them. For 
each pixel recognized as scratched, the suppression 
algorithm defines a 3 x 3 neighborhood around the 
defective pixel. It searches that neighborhood for 

35 any non-defective pixels, and it keeps enlarging the 
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neighborhood until it finds at least one* 
non-defective pixel in the neighborhood. It then 
computes the median of all the non-defective pixels 
in the neighborhood, and replaces the original 
5 defective pixel by this median • This algorithm is 
repeated for each defective pixel. 

The particular algorithm used for scratch 
suppression is not so important as the fact that 
scratches may be detected. 

10 PREFERRED SYSTEM OF THE INVENTION (FIG. 7) 

In accordance with the foregoing detailed 
description, the complete preferred system of the 
invention is illustrated in Fig. 7. A scene viewed 
through camera optics is recorded on film 100 as an 

15 image 102. In addition, several calibration patches 
104 are recorded on another portion of the film 100 
using calibration light source 105 and a gray scale 
step tablet 110. As described in more detail in 
connection with Fig. 1 previously herein, joint 

20 positive and negative developers 120, 122, 124 

generates a joint positive and negative color image 
from the film 100 corresponding to the recorded 
scene image 102 and the calibration patches 104. 
The first step in the process of the 

25 invention is for digital scanner 200 to scan the 
plural calibration patches 104. As described 
previously, each patch is first scanned through a 
red filter and then through a green filter so that 
it separately scans negative and positive images, 

30 An isolation switch 400 routes the digital scanner 
data of the calibration patches to a histogram 
calibration patch processor 402 which forms plural 
histograms in accordance with the previous 
description given in connection with Fig. 2. 

35 Specifically, if there are six calibration patches 



wo 90/01240 



PCT/US89/030S2 



-25- 

within the calibration patch area 104 on the film 
100/ then the histogram processor 402 produces six 
histograms or probability distribution "clouds". An 
interpolation processor 404 then interpolates the 
5 six clouds to form 256 probability distribution 
clouds corresponding to 256 exposure values lying 
within the latitude of the film 100, in accordance 
with the process described previously in detail in 
connection with Fig. 5. These 256 probability 

10 distribution clouds are each designated P(Dp, 
D^|Ej), where the subscript ^ runs from 1 
through 256. 

Next, a large archive of old photographs 
500 representative of the type of scene which is 

15 expected to be recorded on the film 100 is scanned 
by a digital scanner (such as the digital scanner 
200) • This scanning results in a histogram of the 
number of times a particular exposure value^ is 
observed. This produces a probability distribution 

20 P(E) that a particular exposure value will be 

received. A Bayes theorem processor 600 performs 
the computation of equation (1) stated above on page 
17 using all of the 256 probability distributions 
obtained from the- interpolation processor 404 and 

25 the probability distribution obtained from the 

histogram processor 504. A best estimate (mean) of 
exposure processor 602 performs the computation of 
equation (2) stated above on page 18 using the 
results obtained from the Bayes theorem processor 

30 600 to provide a best estimate of exposure for each 
joint value of positive and negative dye density 
Dp, D^. These results are stored in a look-up 
table 300 in accordance with the detailed 
description given previously in connection with Fig. 

35 4. 
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The system is now prepared to process the 
image 102 on the film 100. The isolation switch 400 
now routes the user image data from the digital 
scanner 200 directly to the look-up table 300. For 
5 each pixel in the image 102, the digital scanner 200 
transmits a joint positive/negative dye density 
Dp, through the isolation switch 400 to the 
address input of the look-up table 300. The look-up 
table 300 matches each joint dye density 0^, 
10 thus received to the corresponding estimated 

exposure value E. The exposure value E thus 
determined is transmitted to a printer 900, which 
generates a corresponding image from the succession 

of likeliest exposure values E thus obtained, 

15 In accordance with the description of the 

scratch suppression algorithm given on page 23 of 

the detailed description, a scratch detection 

processor 700 determines from the probability 

distributions generated by the Bayes Theorem 

20 processor 600 which joint dye densities D , D 

p n . 

correspond to unacceptably low probability 
amplitudes for all possible exposure values E. A 
scratch detection. processor 700 causes all such 
; /unacceptable joint dye density values to be flagged 
25 in the look-up table 300. Thus, whenever such an 
unacceptable joint dye density value is received by 
the look-up table 300 from the digital scanner 200 
(through the isolation switch 400), a flag is 
automatically raised so that the corresponding 
30 exposure values stored in the look-up table is not 
transmitted, but instead a special "flag" value is 
transmitted, and a scratch suppression processor 702 
is activated. The missing exposure value is 
provided by the scratch suppression processor 702 as 



35 
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follows: For each pixel recognized as scratched, 
the suppression processor 702 defines a 3 x 3 
neighborhood around the defective pixel. It 
searches that neighborhood for any non-defective 
5 pixels, and keeps enlarging the neighborhood until 
it finds at least one non-defective pixel in the 
neighborhood. It then computes the median of all 
non-defective pixels in the neighborhood, and 
replaces the original defective pixel by this 

10 median. The scratch suppression processor 702 

operates in the same manner for each flag recognized 
from the look-up table 300. 

In summary, the scratch detection processor 
700 examines all the probabilities being computed in 

15 Bayes theorem processor 600, and takes note of any 
probabilities below threshold ("unreasonable- 
values). It then replaces the likeliest exposure 
value in the look-up table 300 with a special value, 
called the "scratch flag" value. Now the look-up 

20 table 300 is ready, and the image is passed from the 
isolation switch 400 to the scratch suppression 
processor 702 through the look-up table 300. The 
scratch suppression processor 702 monitors the 
values as they pass through. Any value which is not 

25 the "scratch flag" value is passed unchanged. When 
a "scratch flag" value is encountered, the scratch 
suppression algorithm is triggered, which grows a 
neighborhood and computes the median, and the 
"scratch flag" value is replaced by the median. 

30 Thus, for scratched pixels the look-up table 300 
does not actually contain "likeliest exposure" 
values, but instead contains a "scratch flag" value. 

Other algorithms 800 may be performed on 
the data before it is transmitted to the printer 

■55 900. Such other algorithms may include various 
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image processes which are well-known in the art for 
enhancing images. 

EXPERIMENTAL RESULTS (FIGS. 8-10) 
The invention has been applied to a real 
5 film. This film is an ordinary, single layer 

emulsion. A glimpse at the experimental results is 
useful. 

Fig, 8 shows the result of histogramming 
the calibration patches for the film, specifically 

10 the probability distributions P(Dp, D^|E) 

obtained from histogramming the calibration patches 
on a real film. The lowest exposure patch is 
plotted in red, the next higher exposure is plotted 
in green, next in blue, next in red, .... Exposures 

15 are separated by 0.3 log exposure (standard 
granularity neutral density step tablet). 
Overlapping locations are plotted in combined color 
(yellow where red and green overlap, cyan where 
green and blue overlap). P(Dp, I^j^IE) for a 

20 single E is named a "cloud," since the edges are 
fuzzy because of the decaying probability 
distribution. Each of the four frames represents a 
slice through the P(Dp, D^^lE) at the percentage 
of peak height shown. The horizontal axis is D^, 

25 the vertical axis is D^. 

The gross features are as expected, but 
many of the details, such as the locus traced by 
successive exposures, and the wide aspect ratio of 
the clouds even at the very low and very high 

30 exposures, are surprising. Fig. 9 shows the result 
of interpolating new clouds (multivariate Gaussian 
distributions). The figure shows 15 clouds for 
clarity. The algorithm was actually run with 256 
interpolated clouds, but this figure shows the case 

35 of only 15 clouds for greater clarity. 
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Fig. 10 illustrates the eight clouds of the 
15 interpolated clouds that correspond to the 
original measured eight clouds. These clouds are 
included in the bottom plot/ but are separated out 
5 in the center plot for easier comparison with the 
top plot- Note that for the frame tagged "> 0.0%", 
the plot simply shows the cutoff for calculating the 
multivariate normal P(Dp, ^^'^^ where the 
probability is very low. In fact, a multivariate 

10 normal distribution is everywhere larger than 0, so 
in principle each cloud fills the entire frame. 

Note that while the choice of using 
multivariate normal distributions to represent the 
clouds is arbitrary, it appears to be a reasonable 

15 fit/ and the results should be fairly insensitive to 
the accuracy of this choice. 

While the invention has been described in 
detail with specific reference to preferred 
embodiments thereof, it is understood that 

20 variations and modifications thereof may be made 
without departing from the spirit and scope of the 
invention. In particular/ there are other ways of 

computing i from D^, D^, but equations 1 and 2 
are the best possible. 

2^ . SOFTWARE IMPLEMENTATION 

A portion of the disclosure of this patent 
document contains computer programs to which a claim 
of copyright protection is made. The copyright 
owner has no objection to the facsimile reproduction 

■^^ by anyone of the patent document or the patent 

disclosure as it appears in the Patent and Trademark 
Office patent file or records, but reserves all 
other rights whatsoever. 

A Fortran program for implementing in a 

35 general purpose computer the noise reduction process 
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described with reference to Figs. 5 and 7 is listed 
in Table I below with reference to the attached 
appendices • 

TABLE I 

5 J I I 



10 



SOFTWARE PROGRAM 


ATTACHED 
APPENDIX 


CORRESPONDING 
FIGURE 


MAKETABLE 
EXPOS 2D 


APPENDIX A 


FIG. 5 


SCRATCH 
SUPPRESS 


APPENDIX B 


FIG, 7 
(700, 702) 
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CLAIMS: 

1. A method for obtaining an enhanced 
image from a photographic film susceptible of 
photographically recording in a frame thereof an 
5 image comprising individual, image pixels 

characterized by individual latent and non-latent 
photosensitive grain counts Dp and 
respectively, said method comprising: 

obtaining a plurality of K distributions of 
10 latent and non-latent grain count pairs D^, 
for K predetermined exposure values, said K 
distributions being at least approximately 
characteristic of said film; 

determining from said K distributions a 
15 likeliest exposure value for each one of a 
plurality of grain count pairs D^, so as to 
generate a corresponding plurality of likeliest 
exposure values; and 

measuring a grain count pair D^, in 
20 a discrete pixel in said film frame and selecting 
the corresponding one of said likeliest exposure 
values so as to generate a set of likeliest exposure 
values representative of the image recorded in said 
film frame. 

25 2. The method of Claim 1 wherein said 

obtaining step comprises: 

exposing a film strip to said K 
predetermined exposure levels so as to form K 
calibration patches in said film strip; 

30 developing the latent and non-latent grains 

in said calibration patches in such a manner that a 
latent grain count and a non-latent grain count 

in discrete pixels therein may be measured 
separately; 

35 measuring latent and non-latent grain count 
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pairs D_ in each discrete pixel in said K 
P n 

calibration patches so as to generate said K 
distributions * 

3, The method of Claim 1 further 

5 comprising the step of printing an image based on 
said set of exposure values representative of said 
image recorded in said film frame, 

4, The method of Claim 1 wherein said 
measuring step is preceeded by the step of 

10 developing the latent and non-latent grains in said 

film frame iii such a manner that a latent grain 

count and a non-latent grain count D in 
p n 

discrete pixels therein may be measured separately. 

5, The method of Claim 2 wherein said film 
15 strip and said film frame are on the same roll of 

film whereby said developing steps are carried out 
simultaneously in said calibration patches and in 
said film frame, whereby latent and non-latent grain 
counts are separately measureable in said film frame 
. 20 as well as in said calibration patches, and whereby 
said K distributions are exactly characteristic of 
: said film. 

6, The method of Claim 1 wherein said 
determining step comprises: 

.- 25 . defining a set of N exposure values E ^ . 

wherein N > K, said N exposure values at least 
partially coinciding with a range of exposure values 
spanned by said K exposure values; 

computing a plurality of N distributions of 
30. grain count pairs D^, for said N exposure 
values from said K distributions; and 

determining from said N distributions said 
likeliest exposure value for each one of said 
plurality of grain count pairs D^, so as to 
35 generate said corresponding plurality of likeliest 
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7, The method o£ Claim 6 wherein said 
computing step comprises interpolating between 
successive ones of said K distributions. 
5 8. The method of Claim 6 wherein said N 

exposure values span an exposure latitude of said 
film strip. 

9 . The method of Claim 1 further 
comprising the step of weighting said distributions 

10 in accordance with a predetermined probability 
distribution of exposure values prior to said 
determining step. 

10. The method of Claim 9 wherein said 
predetermined probability distribution of exposure 

.15 values defines the probability of any pixel in said 
frame being exposed to a given exposure value. 

11. The method of Claim 1 further 
comprising the step of detecting for each pixel of 
said frame whether a scratch is present after said 

20 second measuring step, said method further 

comprising the step of foregoing said selecting step 
for any pixel in which said detecting step step 
detects a scratch and constructing an alternative 
exposure value by interpolating between the exposure 

25 values of neighboring pixels. 

12. The method of Claim 11 wherein said 

detecting step comprises determining the probability 

of occurence of the grain count pair D^, D_ 

^ n p 

measured in each pixel of said frame and flagging 
30 the existence of a scratch whenever said probability 
is below a predetermined threshold value. 

13. The method of Claim 12 wherein said 
probability of occurence is computed for each one of 
a set of exposure values spanning an exposure 

35 latitude of said film strip. 
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14. The method of Claim 7 wherein said 
interpolating step comprises: 

computing a mean value and a variance of 

and and computing a covariance of D 
n p ^ n 

5 with for each one of said K distributions so as 
P 

to define K magnitudes of said mean value, variance 
and covariance corresponding to said K exposure 
values ; 

interpolating between said K magnitudes at 
10 said N exposure values so as to define interpolated 
versions of each one of said mean value, variance 
and covariance at said N exposure values; and 

constructing an interpolated Gaussian 
probability distribution of grain count pairs D^, 
15 for each one of said N exposure values from the 
corresponding interpolated versions of said means 
values, variances and covariance for each respective 
one of said N exposure values. 

15. The method of Claim 8 wherein each of 
20 said K exposure values is defined by an m-bit binary 

word and wherein N = 2™, whereby to utilize the 
full range and resolution with which said exposure 
values are defined. 

16. The- method of Claim 9 wherein said 

25 weighting step computes a normalized weighted sum of 
said K distributions with said predetermined 
probability distribution in accordance with Bayes 
theorem. 

17. The method of Claim 1 wherein said 

30 developing step comprises developing said latent and 
non-latent grains sequentially using a different 
color dye for each. 

18. The method of Claim 17 wherein said 
latent and non-latent grains are developed with 

35 positive and negative dyes, respectively, or 
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negative and positive dyes, respectively. 

19, The method of Claim 6 further 
comprising the step of weighting said distributions 
in accordance with a predetermined probability 

5 distribution of exposure values prior to said 
determining step. 

20, The method of Claim 19 wherein said 
predetermined probability distribution of exposure 
values defines the probability of any pixel in said 

10 frame being exposed to a given exposure value. 

21, A system for obtaining enhanced images 
from photographic film in which there has been 
exposed a film frame to a scene so as to record an 
image therein, said system comprising: 

15 means for developing the latent and 

non-latent grains in said film frame in such a 
manner that a latent grain count and a 
non-latent grain count D^^ in discrete pixels 
therein may be measured separately; 

20 means for storing a plurality of K 

distributions of latent and non-latent grain count 
pairs Dp, for K predetermined exposure 
values, said K distributions being at least 
approximately characteristic of said film; 

25 means for determining from said K 

distributions a likeliest exposure value E^^ for 
each one of a plurality of grain count pairs D^, 

so as to generate a corresponding plurality of 
likeliest exposure values; and 

30 means for measuring a grain count pair 

Dp, D^ in each discrete pixel in said film frame 
and selecting the corresponding one of said 
likeliest exposure values so as to generate a set of 
likeliest exposure values representative of the 

35 image recorded in said film frame. 
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22. The system of Claim 21 further 
comprising: 

means for developing latent and non-latent 
grains in K calibration patches in a strip of film 
5 which have been exposed to said K predetermined 

exposure levels in such a manner that a latent grain 
count Dp and a non-latent grain count in 
discrete pixels in said calibration patches may be 
measured separately; 

10 means for measuring latent and non-latent 

grain count pairs D^, in each discrete pixel 
in said K calibration patches so as to generate said 
K distributions, and for transmitting said K 
distributions to said means for storing. 

15 23. The system of Claim 22 further 

comprising means for developing latent and 

non-latent grains in said film frame in such a 

manner that a latent grain count D and a 

P 

non-latent grain count in discrete pixels 

20 therein may be measured separately. 

24. The method of Claim 22 wherein said 
film strip and said film frame are on the same roll 
of film whereby said K distributions exactly 
characterize said^ film. 

25 25. The system of Claim 21 further 

comprising means for printing an image based on said 
set of exposure values representative of said image. 

26. The system of Claim 21 wherein said 
means for developing said grains in said calibration 

30 patches and said means for developing said grains 
in said film frame are one and the same, 

27. The system of Claim 21 wherein said 
means for determining comprises: 

means for computing a set of N 

35 distributions of grain count pairs D , D for a 

n p 
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plurality of N predetermined exposure values / 
wherein N > K, said N exposure values at least 
partially coinciding with a range of exposure values 
spanned by said K exposure values; and 
5 means for determining from said N 

distributions said likeliest exposure value for each 
one of said plurality of grain count pairs D^, 

so as to generate said corresponding plurality 
of likeliest exposure values. 

10 28, The system of Claim 27 wherein said 

means for computing comprises means for 
interpolating between successive ones of said K 
distributions . 

29* The system of Claim 27 wherein said N 

15 exposure values span an exposure latitude of said 
film. 

30. The system of Claim 21 further 
comprising means for weighting said distributions in 
accordance with a predetermined probability 

20 distribution of exposure values, said means for 
weighting being connected to said means for 
determining . 

31. The system of Claim 30 wherein said, 
predetermined probability distribution of exposure 

25 values defines the probability of any pixel in said 
frame being exposed to a given exposure value, 

32. The system of Claim 21 further 
comprising means for detecting for each pixel in 
said frame whether a scratch is present / said means 

30 for detecting being connected to said means for 

measuring, said system further comprising means for 
foregoing said selecting step for any pixel in which 
said detecting step step detects a scratch and 
constructing an alternative exposure value by 

35 interpolating between the exposure values of 
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33. The system of Claim 32 wherein said 
means for detecting comprises means for determining 
the probability of occurence of a respective grain 

5 count pair D^, measured in a respective pixel 
in said frame and flagging the existence of a 
scratch whenever said probability is below a 
predetermined threshold value. 

34. The system of Claim 33 wherein said 
10 probability of occurence is computed by said means 

for detecting for each one of a set of exposure 
values spanning an exposure latitude of said film. 

35. . The -system of Claim 28 wherein said 
means for interpolating comprises: 

15 means for computing a mean value and a 

variance of and and computing a covariance 
of with Dp for each one of said K 
distributions so as to define K magnitudes for each 
mean value, variance and covariance corresponding to 

20 said K exposure values; 

means for interpolating between said K 
magnitudes at said N exposure values so as to define 
interpolated magnitudes of said mean value, variance 
* and covariance at* said N exposure values; and 
; 25 means for constructing an interpolated- 

Gaussian probability distribution of grain count 
pairs Dp, D^ for each one of said N exposure 
values from the interpolated magnitudes of said 
means value, variance and covariance for each 

30 respective one of said N exposure values. 

36. The system of Claim 29 wherein each of 
said K exposure values is defined by an m-bit binary 
word and wherein N «= 2"^, whereby to utilize the 
full range and resolution with which said exposure 

35 values are defined. 
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37. The system of Claim 30 wherein said 
means for weighting computes a mormalized weighted 
summation in accordance with Bayes theorem of said K 
distributions with said predetermined probability 
5 distribution, 

38* The system of Claim 21 wherein said 
means for developing comprises means for developing 
said latent and non-latent grains sequentially using 
a different color dye for each, 
10 39. The system of Claim 38 wherein said 

latent and non-latent grains are developed with 
positive and negative dyes, respectively, or 
negative and positive dyes, respectively. 

40. The system of Claim 27 further 
15 comprising means connected to said means for 

determining for weighting said distributions in 
accordance with a predetermined probability 
distribution of exposure values, 

41. The system of Claim 40 wherein said 
20 predetermined probability distribution of exposure 

values defines the probability of any pixel in said 
frame being exposed to a given exposure value. 

42. A system for producing an enhanced 
image from film which has been exposed to a user 

25 image and to a plurality of K predetermined exposure 
levels in K corresponding calibration patches, said 
system comprising: 

latent grain developing means for producing 
a dye site of a first color at each latent grain in 

30 said film; 

non-latent grain developing means for 
producing a dye site of a second color at each 
non-latent grain in said film, whereby said latent 
grain developing means and said non-latent grain 

35 developing means together produce a two-colored 
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developed images- 
scanner means for viewing said two-colored 
developed image as an array of N picture elements 
(pixels) and for determining, for each of said K 
5 calibration patches, a dye amount in each of 

said N pixels of said first color, and a dye amount 
of said second color in said developed image; 

first processor means responsive to said 
scanner means for producing, for an i one of 
10 said K calibration patches/ a distribution ^i^^f 



plE^) of a number of pixels having a particular 



joint dye amount D , D , so as to produce an 
ensemble comprising a plurality of K distributions; 

second processor means responsive to said 
15 first processor means for computing, for plural 
joint dye counts 

'^n' Dp in a range, a likeliest exposure value 

E using a normalized weighted sura of products of 
each of said distributions PCD^^^ ^pl^i^ with a 
corresponding exposure level E^; 

means connected to said second processor 

means for 

storing a look-up table correlating each exposure 

value E with a corresponding joint dye amount D^, 
25 Dp; 

third processor means responsive to said 
scanning means for determining from said look-up 
table a likeliest exposure value corresponding to a 
given joint dye count D^, for each pixel in 
30 said user image so as to accumulate an array of 

exposure values corresponding to an enhanced version 
of said user image. 

43. The system of Claim 42 wherein said 
first processor means comprises: 



35 



wo 90/01240 



PCT/US89/03052 



-41- 

fourth processor means for interpolating 
among K of said distributions ^pl^i) 
where 1 ^ i ^ K so as to produce, as said 
ensemble of distributions, M interpolated 
5 distributions ^(0^' ^p'^j^ where 1 ^ j ^ M 
and M > and 

wherein said second processor means 

computes a normalized weighted sum: 

_ M 
10 E (Dp, Dn) « I Ei P(Dp, Dn | Ei) P(Ei)/ 

iel 

M 

r P(Dp, Dn I Ei) P(Ei), 
i=l 

15 44, The system of Claim 43 wherein said 

fourth processor means comprises: 

means for computing, from said K 
distributions, K corresponding values of: 



20 



Dp, a mean value of D^, 



D^, a mean value of D^, 
Spp, a variance of D^, 



25 S^^, a variance of D^, and 

nn n 

Spj^, a covariance of and D^; • 
means for interpolating M points among: 
K mean values D^, 

K mean values D^, 

K variances , 
nn 

K variances S , 

PP 

K covariances S ; 

pn 

35 means for computing/ from said M 

interpolated points, said M interpolated 
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distributions P(D , D^Ie.)* 

n P 3 

45. The system of Claim 42 further 
comprising scratch suppression means, said scratch 
suppression means comprising: 

scratch detection means responsive whenever 

a pixel is sensed in said user image having a joint 

dye count D^, D_ which, when substituted into 
n p 

probability distributions for said exposure levels, 
produces a number below a predetermined threshold, 
for noting said pixel; 

scratch removal means responsive whenever 
said scratch detection notes a given pixel in said 
user image, for ignoring 

a corresponding likeliest exposure value E stored in 
said means for storing a look-up table, and instead 
generating a likeliest exposure value by obtaining 
from said means for storing a look-up table an 
average exposure value of pixels surrounding said 
given pixel but not including said given pixel. 

46. The system of Claim 42 further 
comprising means for scanning a plurality of known 
images and determining therefrom an a priori 
probability distribution of receiving a given 
exposure level for a range of exposure levels 
spanning an exposure latitude of said film, and 
wherein each probability distribution in said second 
processor means is weighted by said a priori 
probability distribution. 

47. A system for producing an enhanced 
image from film which has been exposed to a user 
image and to a plurality of K predetermined exposure 
levels in K corresponding calibration patterns, said 
system comprising: 

latent grain developing means for producing 
a dye site of a first color at each latent grain in 
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said film; 

non-latent grain developing means for 
producing a dye site of a second color at each 
non-latent grain site in said film, whereby said 
5 latent grain developing means and said non-latent 
grain developing means together produce a 
two-colored developed image; 

scanner means for viewing said two colored 
developed image; 

10 first processor means responsive to said 

scanning means for determining a probability 
distribution ^pl^i^ ^^^^ ^ joint dye 

amount D^, in one of said pixels was produced 
by an exposure level in a range corresponding 

15 to said K predetermined exposure levels so as to 
produce a plurality of probability distributions 

P(D^. DpiE.); 

second processor means responsive to said 

first 

20 processor means for determining a likeliest exposure 

level E for a given joint dye count D^, D^, as a 
normalized weighted sum of products of each of said 
distributions ^pl^i) ^ corresponding 

one of said exposure levels E., for plural dye 
25 counts D^, lying within a range; 

look-up table means connected to said 

second 

processor means for storing a likeliest exposure 

E with a corresponding joint dye count D^^, D^; 

means for addressing said look-up table 
means with a joint dye count D^, Dp from a pixel 
in said user image sensed by said scanner means so 
as to obtain a corresponding likeliest exposure 
value so as to accumulate an array of exposure 

35 
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values corresponding to an enhanced version of said 

user image. 

48. The system of Claim 47 wherein said 

first processor means comprises; 

5 third processor means for interpolating 

among K of said distributions P(D^, D^Ie.) 

n p ' 1 

where 1 S i ^ K so as to produce, as said 
ensemble of distributions, M interpolated 
distributions ^p'^j^ where 1 ^ j S M 

10 and M > K, and 

wherein said second processor means 
computes a normalized weighted sum: 
_ M 

E (Dp, Dn) « Z Ei P(Dp, Dn | Ei)P(Ei)/ 



15 



25 



i-1 
M 

2 P(Dp, Dn I Ei)P(Ei) 



i-1 



49. The system of Claim 48 wherein said 
20 third processor means comprises: 

means for computing, from said K 
distributions, K corresponding values of: 



Dp, a mean value of D^, 

D^, a mean value of D^, 

Spp, a variance of D^, 

S„„, a variance of D , and 
nn n 

^pn' 3 covariance of and D^; 
30 . , means for interpolating M points among; 

K mean values D^, 

K mean values D^, 

K variances S , 

nn 

35 K variances S , 

PP 
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K covariances 3^^^; 
means for computing, from said M 
interpolated points, said M interpolated 

distributions P(Djj/ ^pl^j>- 
5 50. The system of Claim 47 further 

comprising scratch suppression means, said scratch 
suppression means comprising:. 

scratch detection means responsive whenever 
a pixel is sensed in said user image having a joint 
10 dye count D^, which, when substituted into 

probability distributions for said exposure levels, 
produces a number below a predetermined threshold, 
for noting said pixel; 

scratch removal means responsive whenever 
15 said scratch detection notes a given pixel in said 
user image, for 

ignoring a corresponding likeliest exposure value 

E stored in said means for storing a look-up table, 
and instead generating a likeliest exposure value by 
obtaining from said means for storing a look-up 
table an average exposure value of pixels 
surrounding said given pixel but not including said 
given pixel. 

51. The* system of Claim 47 further 

25 comprising means for scanning a plurality of known 
images and determining therefrom an a priori 
probability distribution of receiving a given 
exposure level for a range of exposure levels 
spanning an exposure latitude of said film, and 

30 wherein each probability distribution in said second 
processor means is weighted by said a priori 
probability distribution. 
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